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ABSTRACT 


PHASE-LOCK LOOP CHARACTERISTICS IN THE PRESENCE 
OP NON-WHITE NOISE FOR BAND-PASS TYPES OF MODULATING SPECTRA 

by 

ROBERT R. HENRY , B.S. , M. S . 

Doctor of Science in Electrical Engineering 
New Mexico State University 
Las Cruces, New Mexico, 1973 
Professor Frank F. Carden, Chairman 


Conventional analytical techniques used to determine and opti- 
mize phase-lock loop (PLL) characteristics are most often based on 
a model which is valid only if the intermediate frequency (IF) fil- 
ter bandwidth is large compared to the PLL bandwidth and the phase 
error is small. In this paper an improved model (called the Quasi- 
Linear Model) is developed which takes into account small IF Filter 
bandwidths and the non-linear effects associated with large phase 
errors. Results predicted by the Ouasi-Linear Model for a second- 
order PLL are compared to experimental results to determine the 
accuracy of the model. 

Due to its simplicity sinusoidal modulation is frequently used 
to make measurements on an experimental system. In the Experimental 
System used in this paper a more realistic modulating signal is 
used. The spectrum of the signal has a band-pass shape with 
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variable upper and lower break frequencies, with the rolloff rate 
determined by a fourth-order Chebyshev Polynominal. A rather com- 
plex measurement system utilizing cross-correlation techniques 
was required to determine the output signal- to-noise ratio, due 
to considerable overlap in the signal and noise spectra. 

By comparison of theoretical and experimental results it is 
demonstrated that the Quasi-Linear Model accurately predicts PLL 
characteristics. This is true even for small IF Filter Bandwidths 
and large phase errors where the conventional model is invalid. 

The theoretical and experimental results are used to draw conclu- 
sions concerning threshold, multiplier output variance, phase error 
variance, output signal-to-noise ratio, and signal distortion. The 
relationship between these characteristics and IF Filter Bandwidth, 
modulating signal spectrum, and rms deviation is also determined. 
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Chapter 1 


PHASE-LOCK LOOP THEORETICAL ANALYSIS 
1.1 Introduction 

The overview of a typical Communication System is shown in 
Figure 1. The source signal is transformed by the Modulator/ 
Transmitter into a form suitable for transmission over the corn- 
communication channel. By the time this transformed signal has 
reached the Receiver/Demodulator it has been corrupted by additive 
channel noise. The function of the Receiver /Demodulator is to 
perform an inverse transformation on the noise corrupted trans- 
formed signal and obtain a signal closely resembling the original 
signal. 

This paper is concerned with a system such as in Figure 1, in 
which Frequency Modulation (FM) is used to perform the transforma- 
tion. The Receiver/Demodulator function is performed by a Phase- 
lock loop (PLL) . The PLL has become more widely used as an FM de- 
modulator in the past several years. This is primarily due to the 
superior noise— rejection properties and threshold extension exhib- 
ited by the PLL. These characteristics are particularly useful in 
aerospace FM communication systems where receivers must operate in 
low carrier-to-noise ratio (CNR) environments. 

In the development of the traditional PLL model, the additive 
channel noise is assumed to be ‘'white". This implies that the in- 
termediate frequency (IF) Filter bandwidth must be large compared 



COMMUNICATION CHANNEL 



Figure L. Communication System Overview 
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to the PLL bandwidth. It is also usually assumed that the phase 
error is small so that the PLL is operating in the linear region. 

In this paper an improved model which takes into account small IF 
Filter bandwidths and the non-linear effect is developed. 

1. 2 The Phase-Lock Loop Model 

Figure 2 is a more detailed block diagram of the Receiver/De- 
modulator and Communication Channel of Figure 1. The noise process 
n(t) is considered to have a spectral density which is "flat" over 
the passband of the IF Filter. It is assumed that the noise pro- 
cess has a Gaussian amplitude distribution with zero mean. This 
allows n(t) to be expressed as given in [1] by 

n(t) = JT~ [n^Ct) sin (toot) + n^Ct) cos(toot)]. (1“1) 

n 1 (t) and ^(t) are independent Gaussian processes of zero mean and 
identical spectral densities which are the same as the spectral den- 
sity of n(t) but translated downward in frequency so that they are 
centered about zero frequency. Thus the spectrum of n^t) and of 
n^Ct) are low-pass extending to 3 rad. /sec. 

The signal input to the multiplier is given by 

ST~ A sin 0 (t ) + n(t) = (A sin [toot + ( t) ] 

+ n 1 (t) sin toot + ^(t) cos toot} (1-2) 

where A is the rms voltage of the carrier, too the unmodulated car- 
rier frequency in rad. /sec., and 0^(t) the carrier phase due to the 
modulating signal. The Voltage Controlled Oscillator (VCO) output 



*/2Asin0(t) 


Noise 

Process 



Figure 2. Rf ceiver /Demodulator and Channel Block Diagram. 
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signal is 

J2~ K cos 0 1 ( t) = JT~ K cos [d) 0 t + 02(t)J (1“3) 

X 

where is the rms voltage of the VCO output, and G^(t) the phase 
of the VCO output. 

From Equations (1-2) and (1-3) it is straightforward to veri- 
fy that the multiplier output is 

x(t) = AK^ sin[0i(t) - 02(t)] “ '^l n l^ t ^ S ^ n 
+ K^n^(t) cos 0 2 ( t) + AK^ sin[2o>ot + 0 1 (t) + 0 2 C t) 3 ( 1_4 ) 

+ K_n. (t) sin[2(jJ 0 t + 02(t)] + K n (t) cos[2u> 0 t + 02<t)]. 

11 XL 

By neglecting the double frequency terms in x(t) , Equation (1-4) and 
Figure 2 yield the following for the VCO input 
t 

e( t ) * K J* [A sin <j>(u) - ^(u) sin 0 2 (u) + n 2 (u) cos 0 2 (u)] 

f (t-u)du U“5) 

where f ^ l cue impulse ~ * ’ * 

({>(t) = 0i(t) - 0 2 (t). (1”6) 

Here the PLL is assumed to be locked at t=Q so that e(0) — 0. Let 
K 2 be the VCO constant (rad. /sec. /volt) , then 

iigli i = K,e(t) + u) 0 . (1-7) 

dt *■ 

Also define a new constant 


K • 


(1-8) 



6 

Now Equations (1-5) through (1-8) yield the phase model of the PLL 
shown in Figure 3. 

Now consider the case in which 02 (t) is small (less than 30°). 
This would be true for a small modulation index and high CNR. Since 

1*2 (t) cos 0 2 ( t ) - n^(t) sin 02 (t) — ► ^(t) A n(t) (1-9) 

as 02 (t) — ► 0, the phase model of Figure 3 becomes that of Figure 4. 
The spectrum of n(t) is low-pass with cutoff frequency $ and spectral 
density the same as the noise process as discussed in the noise rep- 
resentation given by Equation (1-1). Viterbi [1] shows that the 
model of Figure 4 is valid for the case in which the IF filter band- 
width is wide compared to the PLL bandwidth. Thus the model given 
by Figure 4 is valid for small phase errors and/or wide IF Filter 
bandwidth compared to PLL bandwidth. 

Since the multiplier output (neglecting double frequency t^rms) 
is 

x(t) = K. [A sin <t>(t) + n(t)l. (1-10) 

one finds that 

x(t) = K x y (t) . (1-11) 

One of the difficulties in analyzing the PY1L represented in 
Figure 4 is that a linear analysis is invalid t tot —.til c(~' 

since sin <£(t) = <j>(t) only for <$>(t) < 30°. ..^i. .. c' 

difficulty is to replace the non-linearity by su v aie..c 

gain which depends on 4>(t). This quasi-linear at .-..b vas first 
applied to the PLL by Develet [2]. Replace the n<'c-linerrity bv ^ 


















9 


linear gain G so that the expected value of the square of the dif- 
ference E {[A sin <KO - G cj>(t)] 2 } is minimized. Assuming that 4»(t) 
is a zero-mean Gaussian process with a 2 = E[(j> 2 (t)] it can be shown 
[3] that 

G = A e" a2/2 . (1-12) 

The assumption that 4>(t) is Gaussian is justifiable only when the 
PLL is operating in the linear region. For large (>30°) <f) the Gaus- 
sian assumption may no longer be valid, but the above procedure 
should be more accurate than replacing A sin 4> (t) with A <f>(t). The 
phase model now becomes that illustrated in Figure 5. The transfer 
function D(s) represents a post— detection filter. 

1.3 Steady State Analysis 

Assuming that the modulation and additive noise processes are 
zero mean and wide sense stationary, it is straightforward to show 
[1] that 4>(t), y(t) , o(t) , 0 i ( t ) , and 0 2 (t) are zero mean and that 
as t (steady state) the covariances are independent of t. Thus 

the corresponding spectral densities are well defined. Since the 
non-linearity is replaced by an equivalent linear gain G, the system 
in Figure 5 is linear and the superposition principle holds. This 
coupled with the hypothesis that the modulation process and the ad- 
ditive noise process are independent allows one to compute the ef- 
fect of the modulation and of the noise independently and then com- 
bine the two. 

Let 0i Go), 02 (oj) , ?(u>), N(<o), 7(e)), and 0(w) represent the 
spectral densities of 0i(t), 02(t), <H C ). n (h)» y(t)j and o(t) 



n(t) 



o 
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respectively. Throughout this paper it is understood that s is 
written in place of j«. It can be shown [4] that the frequency 
response of the PLL is given by 


H(s) A 


02 (s) _ GK F(s) 
0i(s) s + GK F(s) 


(1-13) 


where 0 2 (s) and 0i(s) are the Fourier Transforms of 0 2 (t) and 0i(t) 
respectively. 

Equation (1-13) and Figure 5 give (see Appendix I) 

$(u>) = |l - H(s) | (to) + |h(s) | N(o)), (1-14) 

G 

Y(w) = G^|l - H(s)|^0i(w) + | 1 - H(s)| N(w), (1-15) 


and 

0(w) = oi 2 |H(s)| 2 D( U )0i(ai) + — |h(s) | D(o))N(u)) . (1-16) 

G 

From Equation (1-14) the variance of the phase error is 

/*•»£ 

T — DO 

= /|l - H(s)| 2 0i(td) \ / |H(s)| 2 N(w) || . (1-17) 

-CO G -CO 

From Equation (1-11) and (1-15) the phase variance of the multiplier 


output is 


= (GK ^) 2 


00 



— QO 


/ |l - HCs)! 2 ^!^) |£ + K 2 / |l - H(s)| 2 N(oj) 


dco 
2 tt * 

(1-18) 


Likewise Equation (1-16) gives the output variance as 
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£ r — , N da> 

= / 0(M) 2? 

— a 

dw 


/ uT|h(s) | 1 D(a>) 0i (oo) |£ + -y / w 2 |H(s)| 2 D(a>) N(a>) 
-oa G —00 


da) 

2tt 

(1-19) 


1.4 Rela tionship Between g 2 and o' 2 

In much of the literature the phase variance is often deter- 
mined by measuring O 2 and multiplying by the appropriate phase de- 


tector constant. From Equations (1—17 » 1—18) it can be seen that 

■ < GK r 


2 / v 2 2. \ __ 

this is true only for the no noise case (here a, « (GK ) 0 ). For 


the general case with noise 

2 


du> 

277 


a 2 - a? = [(GK,) 2 - 1] r 1 1 - H(s)| 2 0i (a>) 

x 9 1 J 

— 00 
oo 

+ / [k 2 |i - H(s )| 2 - -||h(s)| 2 ]n<w) . 

—00 L ’ G 

Using the fact that 


(1-20) 


|1 - H(s) | 2 = |H(s) | 2 


or 


(GK) 2 | F(s) 


( 1 - 21 ) 


Equation (1-20) becomes 


- a 2 = [(gk x ) 2 - 1] f |l - H(s) | 2 0i (w) 


00 

4 f [-4^ j " A jHCs) I 2 NCto) 

G 2 1 Lk 2 |F(s)| 2 J 


— n da) 
2tt 


( 1 - 22 ) 


Notice that the relationship involves the PLL parameters as well as 


the noise spectral density. 



Chapter 2 


THE EXPERIMENTAL SYSTEM 


2.1 System Overview 

An overview of the Experimental System is shown in Figure 6. 
Both the noise generator used to produce the signal and the one 
used to produce the additive channel noise have an essentially 
"flat” spectral density over the bandwidth of the corresponding 
filters which follow them (B(s) and I(s) respectively). The cen- 
ter frequency of the FM transmitter was chosen to be 455 kHz* so 
that the IF Filter is centered at 455 kHz. The noncommercial equip- 
ment schematics are shown in Figures 7 and 8. 


2.2 The Band-pass Filter 

The Band-pass Filter used to shape the "white” noise is the 
General Radio Type 1952 Universal Filter. The cutoff frequencies 
are selectable with the response being that of a 4-pole Chebyshev 
filter. Such a response for the low-pass section is given by 


|G(s)| 2 = 


1 + e 2 T 4 


ft) 


(2-1) 


where T is a 4th order Chebyshev Folynominal and the peak- to peak 
4 


ripple is given by 


1 - 


r — r 

1 + e 2 


( 2 - 2 ) 


The parameters as determined from comparison of the measured 
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Figure 6. Experimental System Overview. 
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Figure 7. Multiplier and Summer Schematics 
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frequency response and theoretical response (see Figure 9) are: 


e 2 = 0.07, f x = 1.2 f 1>3dB , f 2 = 0.895 f 2j3dB (2-3) 

where f n is the selected cutoff frequency for the high-pass fil- 

ter section, and f the corresponding theoretical break frequency. 

A similar statement applies to f^ ^2 ^ or t * ie low-pass section. 

By cascading the low-pass and high-pass sections the desired trans- 
fer function for the band-pass filter is obtained. The result is 
(see Appendix II) 

4 

B(s) = q — “ — * (2-b) 

ggS + g ? s + + g ± s + g 0 

The constants g Q through g g are functions of the upper and lower 
break frequencies and are given by Equations (II-5) and II-9). 


2.3 The FM Transmitter 

Consider the effect of applying a sinusoidal signal of fre- 
quency f Hz and zero-to-peak voltage of V to the input of the FM 
transmitter. Then the zero-to-peak deviation (Hz) Af of the carrier 
is 


2TtAf = XV. 


(2-5) 


It is straightforward to verify that the zero-to-peak phase devia- 
tion is 


From Equation 


X = 


3 = Af/f 

m 

(2-5) and (2-6) 

2tt 3 f 

m 

V 


= Ao)/(jJ 

m 


rad/sec 

volt 


(2-6) 


(2-7) 




Frequency(Hz) 

Figure 9. Band-pass Filter Frequency Response. 
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It is well known [5] that for sinusoidal modulation the am- 
plitude of the carrier is proportional to J Q (3) where J Q ( ) 1S 
the Bessel function of the first kind of order zero. Thus for 
0 = 2.41 the carrier component vanishes. 

The above 'mentioned modulation was applied to the FM Trans 
mitter and V increased until the carrier component of the Trans- 
mitter vanished as observed on a Spectrum Analyzer. The value of 

V and f were then recorded. Then Equation (2-7) with 3 ” 2.41 was 
m 

used to determine X. Values of f m from 400 Hz to 2000 Hz were used 
with less than 0.1% difference in the corresponding values of X. 

X was determined to be 


X = 2tt(201000) 


rad / sec 
volt 


( 2 - 8 ) 


2.4 The IE Filter 

The measured frequency response of the four IF Filters used m 
the Experimental System is shown in Figure 10. The Figure is actu- 
ally a plot of l(s) translated down in frequency so that the center 
frequency becomes zero. It is convenient for future calculations 

to define 

C(s) A l{j 27T(f - 455000)} (2-9) 

since it is the translated version of I(s) that is needed. 

The frequency response of B(s)^p determined in Appendix II with 
the appropriate break frequency fits the frequency response of C(s) 
reasonably well as illustrated in Figure 10. Therefore from 
Equation (II-6) 



|l(joj)| (dB) 
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( 2 - 10 ) 


P 4 s 4 + P 3 s + p 2 s + p x s + f o 


P 3 = f 3^' 


P 2 - f 2 /B J 


P 4 - 1/3 


( 2 - 11 ) 


and B is the break frequency in rad/sec. The values of B for the 
four mechanical filters that could be used as the IF filter are 
given in Figure 10 along with the theoretical and measured frequency 

responses . 

2.5 The Post-Detection Filler 

The Post-Detection Filter used has a first order low pass re 
sponse given by 


, . D 

D < s > - b^n 


(2-12) 


where b is the break frequency in rad/sec. The filter is used to 
eliminate output noise outside of the signal spectral characteristic 
and thus enhance the output signal- to-noise-ratio. 


2.6 The Multiplier 

The Multiplier is shown in Figure 7 and is part of the EEP XR- 
S200 Multifunction Integrated Circuit. The resistor-capacitor com- 
bination R -C £ filters out the double frequency terms of the multi- 
plier output. The multiplier output measurements were made at point 
0of Figure 7. With the loop switch set to the "open loop" position 
and with no noise, the multiplier output was observed on an oscil- 
loscope for various offsets in frequency between the VCO and FM 
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transmitter. The observed waveform was sinusoidal with constant 
amplitude for offset frequency differences up to 5 kHz. Since the 
highest modulating frequency used in the system was about 1.5 kHz, 
the multiplier response is essentially flat over the frequencies of 
interest. The zero-to-peak voltage output of the multiplier output 
was measured to be 1.4 volts. Thus from Equation (1-10) with n(t) = 0 

K^A - 1.4(volts)^ (2-13) 

The voltage of A measured at (IT) in Figure 7 was found to be 0.776 
volts rms. Thus from Equation (2-13) 

K, = * 1,8 volts (2-14) 

1 0.776 

2.7 The Loop Filter and PLL Transfer Function 

A very common loop filter used in PLL s is one with an integral 
plus proportional 11 transfer characteristic. This type filter is 
used in the Experimental System and has a transfer function given by 

F(s) - 1 + | . (2-15) 

It follows from Equation (1-13) that 

H(s) - f ^ (2-16) 

s z + GKs + aGK s + x.,s + x Q 

where 

x Q = GKa, x 1 « GK. (2-17) 

Since the largest exponent of s in the denominator of Equation 
(2-16) is 2 the PLL is classified as second order. The actual loop 
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filter is shown in Figure 8. It can be shown [4] that for a large 
operational amplifier gain this active filter has a transfer function 

given by 


x _i + ^o_ 

T. T 8 


ST ^ + 1 
ST 


2 + 1 T 2 r 1 1 

i = — 1 + 

1 T lL T 2 S . 


(2-18) 


where = R-^C, ~ ^2^ * ~ 

Comparing Equations (2-15) and (2-18) one finds that 

a = 1 /t 2 ( 2 ~ 20 > 

From Equations (1-13) , (2-15) and (2-18) it can be seen that the 
loop gain K is modified by the factor T^. Gardner [4] shows 
that 

"n = ^ and C - \ « n . C 2 ’ 21 ) 

where 0) is the natural frequency and £ the damping factor for the 
n 

second order system. In this system the following values for Ui^ 
and C were chosen 

03 = 271(860), £ = 0.707 (2-22) 

n 

The values of R , R 2 , and C were chosen to give the approximate 
values for d> n and £. Then the loop gain was changed (by varying R g ) 
until the measured frequency response agreed closely with the the- 
oretical response [4], 

The Post-Detection filter was chosen so that 

b = 2 tt(1300) . (2-23) 

The overall measured frequency response of the PLL is shown in Figure 
11. Also included for comparison is the theoretical response for 




Frequency (Hz) 


Figure 11. PLL - Post-Detection Filter Frequency Response. 


NJ 
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the parameters given in Equations (2 22) and (2-23). 

2.8 CNR, SNR, and Voltage Measurements on the System 

CNR measurements were made with a rms voltmeter at point ® in 

Figure 7. This is equivalent to measuring r(t) in Figure 6. The 

carrier rms voltage (V ) was measured with no modulation applied to 

the FM transmitter and no channel noise added. A calibration chart 

was then made for each IF Filter which gives the voltage at point® 

(V ) for a given voltage at the channel noise generator output with 
n 

the carrier turned off. From this calibration chart the CNR was 
determined by reading the rms voltage at the channel noise gener- 
ator output. That is 

CNR = 20 log 1Q ( v c / v n > (2-24) 

where V was determined from the calibration chart and the noise 
n 

generator output reading. 

The multiplier output voltage was measured with a rms voltmeter 
at point® of Figure 7. This point corresponds to the measurement 
of x(t) in Figure 6. Thus the rms voltage reading squared gives 

the multiplier output variance. 

The measurement of the output SNR presented a problem that is 
not encountered when the modulating spectrum is narrow compared to 
the loop bandwidth. The solution required the use of a relatively 
complex and powerful measurement device, namely a HP 2116 computer 
connected to two analog-to-digital converters via a HP 2115 computer. 
Quite often SNR measurements are made by using filters to filter 
out" most of the noise or the signal to determine the 
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signal - to - noise ratio. The fact that the signal and noise output 
spectrum overlap considerably precludes the use of such a technique 
in this system. The inherent non-linearity of FM and the PLL also 
preclude the use of superposition to measure each effect separately. 

Since frequency-domain diversity characteristics could not be 
exploited, time-domain differences were utilized. The cross cor- 
relation (R (r)) between the input signal (s(t)) to the FH Trans- 
so 

mitter and the output signal (o(t)) was the fundamental characteris- 
tic used to determine the SNR. Figure 12 is a block diagram of the 
measurement system. The pulse generator was adjusted so that both 
s ( t ) and o ( t ) were sampled at a 25 kHz rate. Since the highest fre- 
quency content of the signal was about 1400 kHz, the sampling rate 
was approximately 18 times the highest frequency. The analog-to- 
digital converters (ADC) were adjusted to quantize the signals into 
a 10 bit word (1024 levels). Since the ADC's required a 0-10 volt 
input, the signals were offset by approximately 5 volts and ampli- 
fied to cover a dynamic range from about 2 to 8 volts. One thou- 
sand samples were transferred to the 2116 computer by first setting 
parameters on the experimental system and then running the computer 
program. 

A block diagram of the computer program used to calculate the 
SNR is shown in Figure 13. Due to the IF Filter and PLL character- 
istics the output signal is delayed with respect to the input signal. 
The approximate delay was determined by observation of the signals 
on an oscilloscope. Then two numbers (ID1, and 1D2) such that 
ID1 < // of samples of delay < ID2 



PULSE GENERATOR 


(one pulse 
every 40 ys) 



Figure 12. Measurement System Block Diagram 









START 


DO 20 
G = G1,G2 


N ~ 1000 


COMPUTE, STORE MSE(G) = 
N-IDM 

F [ o(t i-IDM ) ' Gs(t i ) ^ 
i=XDM+l 


INPUT: ID1, ID2 , 
G1,G2 


20 


READ, STORE ADC SAMPLES 


DO 10 

ID = IDl,ID2 


FIND GMIN SUCH 
THAT MSE(G) IS 
MINIMUM FOR 
G = GMIN 


COMPUTE .STORE R (ID) = 
7 so 

N-ID 

N ' s(t.)o(t. TT J 

N-ID /j V x 7 l-ID 

i=ID+l 


N-IDM 

S 2 2 
GMIN s (t i ) 

[+1 


FIND IDM SUCH THAT 
R so (ID) IS MAXIMUM 
FOR ID =IDM 


SNR =10L0G[SIG PWR/MSE(GMIN)] 


Figure 13. Measurement System Computer Program Flow Chart, 
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were read into the computer. The program then determined the delay 

(IDM) for which the cross correlation R (t) was maximum, in the 

domain ID1 to ID2. Next the mean-square-error (MSE(G)) between 

s(t) and o(t) was determined for the input signal delayed by IDM 

samples and the output signal multiplied by a gain G, for G1 _< G 

< G2. The values G1 and G2 were previously input into the computer 

such that a plot of MSE(G) showed a distinct minimum. The output 

SNR was then calculated as shown in Figure 13. In summary then, the 

SNR was calculated by delaying the input signal s(t) by IDM samples 

such that the cross correlation R (t) was maximum, and multiplying 

so 

the Input signal by a constant G such that the mean-square-error 
MSE(G) was minimum. 

2. 9 Experimental System Measurement Errors 

The ideal Experimental System would have components which are 
described exactly by their mathematical representation. However, 
this is not usually the case. For example, a resistor may not exhibit 
a linear relationship between voltage and current although it is 
modeled this way. Thus, there are errors introduced by the Experi- 
mental System which are not accounted for by the mathematical model. 

In well designed systems these errors are kept small so that they 
will not interfere with the primary purpose of the system. 

The Measurement System illustrated in Figure 12 introduces some 
error in the SNR measurement. This is due to the slight differences 
between and the quantization errors of the two Analog-to-digital Con- 
verters, as well as the finite sampling rate. To determine the 



magnitude of these errors, the same analog signal was applied to the 
inputs of both ADC’s and the output SNR calculation observed. For 
signals in the frequency and amplitude range of those used in this 
paper the SNR reading was greater than 30dB. Ideally the reading 
would be infinite since the two signals are identical and the mean- 
square-error is zero. Thus the "noise" introduced by the measure- 
ment system is at least 30dB below the signal being measured. 

The main source of error in the Experimental System is the 
"phase jitter" introduced by the VCO. With no modulation applied 
to the system the phase error should be zero for the ideal system. 
However, oscillator instability in the VCO and Transmitter cause 
"noise" at the output of the multiplier. The magnitude of these and 
other errors were determined as follows. The Transmitter was modu- 
lated by the signals with spectra as illustrated in Figure 9 using 
a small modulation index so that the PLL was operating in the "linear 
region. The Measurement System was then used to calculate the output 
SNR. This value ranged between 18 and 20dB for the various signals. 

The combined "noise" of the Measurement and Experimental Systems 
is approximately 20dB below the signal being measured. T he majority 
of this noise is caused by "phase jitter" of the VCO and Transmitter. 
Thus SNR measurements below 20dB should be accurate in the sense 
that it is the Experimental System characteristics and not system 
design errors which are dominant in causing the reading. 



Chapter 3 


THEORY APPLIED TO THE EXPERIMENTAL SYSTEM 


3.1 RMS Deviation and Noise Power 

Since the input signal to the FM transmitter is essentially a 
white noise process shaped by a band-pass filter, the spectrum of 
s(t) assumes the filter spectral shape and is given by (see Equation 


(2-4)) 


S (s) = R B(s) = 


Rs 


7 

+ g 7 s + 


(3-1) 


+ gjS + g 0 


where R is a constant which is related to the rms voltage of s(t) 
via „ 


v 2 = R 2 f B(u>) 
rms J 2tt 


where B(w) = |B(s)| 2 (3-2) 


Let X be the FM transmitter constant in rad/ sec. /volt, then 

DO 

(RMS Deviation) 2 = (XR) 2 J B(w) |f ‘ (rad/sec) 2 . (3-3) 

— 00 

The phase of the carrier is given by the integral of the frequency 
of the carrier with the appropriate initial conditions. Assuming 
zero initial conditions, Equation (3-1) gives 


01 (s) = fs(s) = RX^. 


(3-4) 


Thus the spectral density of the carrier phase due to the signal s(t) 
is 


0 (w) 

l 


(RX) 2 B(d)) 

2 

u) 


(3-5) 


1 
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f 

The noise process n(t) has a low-pass spectrum with cutoff 
frequency 0 rad/sec. (see Section 1.2) Since the shape is 
determined by the IF Filter one finds from Equation (2-10) that 


N(s) = TC(s) 


P 4 S + p 3 5 


( 3 - 6 ) 


T is a constant which is related to the rms noise voltage into the 
multiplier via co 


noise power = v^ mg - 2 J N(O)) 2tt 


where 


N(03) = |N(s) 


( 3 - 7 ) 


( 3 - 8 ) 


3.2 Folynominal Form of the Spectra 

The general equations developed in Chapter 1 for the PLL spec- 
tral densities are now specialized for the Experimental System of 
Chapter 2. From Equations (2-16), (2-4), and (3-4) 


[1 - H(s)][Gi(s)] = 


4 


L s + X x s + X Q J 


RXs' 


„ .8 A 
* & 8^ 


.+ g lS 


+ g 0- 


RXs' 


h io sl ° + V 9 + " 


( 3 - 9 ) 


. .+ h^s + 


where 

h 0 = x o g o 
h l = x l g 0 + x o s i 

h i = g i-2 + x l 8 i-l + X 0 S i for 2 < i < 8 
h 9 = g 7 + x l g 8 


h 10 ~ s 8 ' 


( 3 - 10 ) 
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Equations (2-16) and (3-6) yield 
[1 - H(s) ] [N(s) ] = 


s 2 1 

T . I 

2 . 

_s + x^s + XqJ 

J> 4 s 4 - + P x s + f 0 J 


Ts' 


e^s +.«..+ e^s 4 - 


(3-11) 


where 


e 0 x 0 £ 0 

e i = x i f o + x o p i 

e 2 + f 0 + X 1 P 1 + X 0 P 2 

e 3 = P 1 + X 1 P 2 + X 0 P 3 


e ; = 


P 2 + Xl P 3 + x 0 p 4 


e 5 = p 3 + x lP/i 


e 6 = p 4 


(3-12) 


From Equations (2-16) and (3-6) 
H(s) N(s) = 


r x i s + x o 

X 1 

Ls 2 + x^s + X Q _ 

_P^S Pj_ s ^ ^ Q J 

T V + 

Xx 0 


e,s +....+ e.s + e A 

6 X U 


Using Equations (2-16), (2-12), and (3-4) 


(3-13) 


sH(s) D(s) 0i (s) = 


X^S + XqS 


S + X^S + Xq 


b + s 


RXs" 


d 3 s + d Q s 


4- d^s^ + ^i s **" ^0 


l^ggS +. . . « 4 * g^s 4 - 

3 


RXs 


g g s +....+ g s + g Q 


(3-14) 


where 


d 0 = x o b 


d 2 = b + X 1 


d 1 = X x b + x 0 


dg - bx^ . 


(3-15) 
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Further algebraic manipulation of Equation 


RX d s 5 

s H(s) D(s) 0i (s) = — — — Jo 

h u s + h 1Q s 


(3-14) gives 


+ d 0 s 


+....+ h^s 4* 


(3-16) 


where 


II 

o 

\ja 

s o d o 




h l 

g o d i 

+ 

g l d 0 


IE 

\si 

g 0 d 2 

+ 

g l d l 

+ 8 2 d , 

II 

1** 

g i-3 

+ 

g i-2 d 2 + 8 

h 9 = 

g 6 + 

g 7 

d 2 + 

s 8 d l 

K 81 

o 

II 

g 7 + 

g 8 

: d 2 


h ll” 

g 8 ' 





3 < i < 8 


(3-17) 


Equations (3-6), and (3-14) give 


s H(s) D(s) N(s) = 


d^s + d Q s 


s d + d^s^ + d^s + d^ 


p y s +....+ p., s + f n 
4 1 U 


T (d 3 s + d Q s) 

- 7 "T T 11 . - 

e^s +....+ e^s + eQ 


(3-18) 


where 




f o d o 


e 5 = p 2 + p 3 d 2 + p 4 d 1 


e, = 


f 0 d l + p i d o 


e 6 = p 3 + p 4 d 2 


" f 0 d 2 + Pi d i + p o d 


11 


2 0 


e„ = 


f 0 + p l d 2 + p 2 d l + p 3 d 0 


e 7 = p 4 


e 4 = P 1 + P 2 d 2 + p 3 d l + P 4 d 0 


(3-19) 



3.3 Simplification of Results 

The following definitions are now made to simplify notation 
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Recalling that s = jw, 
INTI A 


/ 1 


+....+ h^s + hg 


2 dw 
2tt 


INT2 A 


cx 

/ 


d 3 s 5 + d Q s 4 


Ih^s"^ +,...+ h^s + hg 


2 dw 
2tt 


INT3 A 


a 

/ 


e^s 4“ « • . .4* e^s e Q 


2 dw 
2tt 


INT4 A 


f I *i s + x o 

J le^s^ +**■•+ + 


2 dw 
2 rr 


INT5 A 


[ L d 3 S -1 

J I e^s^ +....+ 


dg® 


e l s + e o 


2 dw 
2 tt 


INT6 A 


/ 1 » , •*, , 

J IggS +*•••+ g-^ + 


2 dw 
2tt 


INT7 A 


a 

/ 


4 j. 

ip 4 s +. 


. .+ p- s + f ( 


2 dw 
2tt 


(3-20) 


By combining Equations (1-17), (3-9), (3-13) and (3-20) one 
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obtains 


c! = (RX) 2 (INT1) + ~ (INT4) 
t nr 


(3-21) 




4>,n 


In the above 0 , is the phase error variance due to the signal while 
4>,s 

2 

a, is the variance due to the noise. 

<j>,n 

In a similar manner Equations (1-18), (3-9), (3-11) and (3-20) 


give 


= (GK X ) 2 (RX) 2 (INTI) + (K^T) 2 (INT3) 


(3-22) 


x, s 


x,n 


while Equations (1-19), (3-14), (3-16), and (3-20) give 

2 2 T 2 

0 = (RX) Z (INT2) + (INT5) 

0 G 2 


o,s 


o ,n 


(3-23) 


2 2 2 2 
where 0 a , and 0 are defined in a manner similar to 

x,s x,n’ o,s o,n 

2 2 

a, and 0, 

(p y S y j n* 


From Equations (3-20) and (3-1) , Equation (3-3) becomes 


(RMS Deviation ) 2 = X 2 R 2 (INT6), 


(3-24) 


while the from Equations (3-7), (3-6), and (3-20) the noise power 
at the input to the multiplier is given by 


noise power = 2T (INT7). 


(3-25) 
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3.4 Computer Implementation of Equations 

A Computer Program was written to evaluate the quantities in 
Equations (3-22) through (3-25). The block diagram of the Program 
is shown in Figure 14, For the input parameters Wi , ( 1)2 (signal 
break frequencies), $(half the IF bandwidth), and DEV (rms voltage 
into FM Transmitter) , the program calculates the indicated values 
for CNR' s from 0 to 27dB in 3dB steps. The output values for each 
CNR is punched on paper tape for later use. 

The value for G was determined by first setting G - A, calculating 
2 2 

a,, then setting G = A exp(-a,/2) as indicated by Equation (1-12). 

2 

Then O. was re— calculated using the new value of G. This process 
was iterated three times to obtain a third-order estimate of G. 

The integrals were evaluated by numerical integration (Simpson's 
Rule) after much difficulty was encountered with analytical tech- 
niques. The equations for the analytical method were quite complex 
and often caused overflow and underflow during computer calculations. 
The integration increment used for the numerical integration was 60 
rad/sec (about 10 Hz). This value was decided on since reducing the 
increment 100 times changed the integrals by less than 0.001%. The 
lower limit of the integration was zero while the upper limit a fre- 
quency which corresponds to the response being 30dB smaller than the 
maximum. The 30dB figure was chosen since setting the value to 90dB 


caused less than 0.01% change in the value. 

The computer used was a HP 2116 programmed using the Fortran IV 
language. Typically about 15 minutes was required to obtain a complete 
























set of output data for a given set of input parameters. Over 95% 
of this time was used to evaluate the integrals by the numerical 


technique . 



Chapter 4 

PLL CHARACTERISTICS IN THE ABSENCE OF NOISE 

4.1 Signal Distortion Due to the IF Filter and ?LL Nonlinearity . 

The spectrum of the transmitter output is a function of the 
signal and deviation of the carrier. For sinusoidal modulation, 
98% of the total power in the FM signal is contained in a band- 
width (B) given by Carson’s rule [5]: 

B = 2(Af + f ) (4-1) 

where Af is the zero-to-peak frequency deviation and f^ is the 
modulating frequency. For a modulating signal with a bandpass 
spectrum the corresponding FM signal bandwidth is difficult to 
determine. Van Trees (pp. 100-104 of [7]) indicates that a 
measure of the bandwidth for a modulating signal with a low— pass 
Gaussian spectrum is 2 a d^., where O is the signal rms voltage 
and d the Transmitter constant. 

From the above it is reasonable to assume that the larger 
the rms deviation the larger the bandwidth occupied by the FM sig- 
nal. Thus increasing the rms deviation and/or decreasing the IF 
Filter Bandwidth causes more of the FM signal to be filtered out 
and results in more distortion in the demodulated output. This 
effect is illustrated in Figures 15 through 18. Note that a 
decrease in IF Bandwidth and/or an increase rms deviation results 
in a lower SNR (larger signal distortion). 

In addition to the signal distortion caused by the IF Filter, 
distortion is also caused by the fact that the PLL is non-linear. 



RMS DEVIATION (HZ) 


Figure 15. SNR as a Function of Deviation for No Noise with f =56Hz,f “440Hz 
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Figure 16. SNR as a Function of Deviation for No Noise with f 82 56 Hz, f 2 « 900 Hz. 
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Figure 17. SNR as a Function of Deviation for No Noise with = 56 Hz , = 1250 Hz. 
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Figure 18. SNR as a Function of Deviation for No Noise with ■ 560 Hz, f ^ 
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This distortion is more pronounced for larger O ^ since the non- 
linearity sin is nearly linear for small but non-linear for 




large (>30°) a^. Now an increase in the ms deviation "stresses" 
the PLL more and results in an increase in 0 Thus as the ms 
deviation is increased the PLL output SNR decreases (distortion 
increases). If the loop "stress" is sufficient a phenomenon known 
as "cycle-slipping" will occur [8,12] resulting in impulse noise at 
the PLL output. A further increase in the ms deviation will nor- 
mally result in a complete loss of phase lock [9]. The "cycle- 
slipping" and "loss of lock" phenomenon was observed as the ms 
deviation was increased for all four signals and four IF Filters 
used in the Experimental System. A summary of distortion due to 
PLL Stress and Finite IF Filter bandwidth is presented in Table 1. 
Data for this Table was obtained from Figures 15 through 18. 

Figures 15 through 18 illustrate that for no IF Filter 
(infinitely wide IF bandwidth) an increase in the signal upper cut- 
off frequency results in a lower output SNR. In obtaining the data 
the deviation was increased until the loop would not maintain lock. 
The ms deviation at which loss of lock occurs (obtained from 
Figures 15-18) as a function of f^ for f^ * 56 Hz Is plotted in 
Figure 19. A non-linear analysis of the second-order PLL for sinus- 
oidal modulation [10,11] has shown that the PLL unlock boundary for 

f < 0.85f is given by: 
m n ° 


Af (f ) = 1.15 f 2 (4-2) 

m n 

where Af is the zero-to-peak deviation, f the modulating frequency, 
and f the PLL natural frequency. Since it is the highest modulating 
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Table 1 

Output Signal Distortion Due to PLL Stress 
and Finite IF Filter Bandwidth for No Noise 


IF FILTER 

signIF^“™ 

(fl - f 2 )^l 

5600 HZ 

3300 HZ 

900 HZ 


1. OdB 

1 . OdB 

1. -8dB ,IF 

56 HZ - 





2. OdB 

2. OdB 

2. -10dB,IF 

440 HZ 

. 




3. OdB 

3. OdB 

3. -15dB ,IF 


1. -3dB,PLL 

1. -5dB,PLL 

1. -lOdB ,IF 

56 HZ - 




900 HZ 

2. -3dB,PLL 

2. -5dB,PLL 

2. -15dB,IF 


3. ~6dB ,IF 

3. -10dB,IF 

3. -18dB, IF 


1. -6dB,PLL 

1. -6dB,PLL 

1. -16dB ,IF 

56 HZ - 





2. -6dB,PLL 

2. -6dB,PLL 

2. -18dB,IF 

1250 HZ 





3. -lOdB.IF 

3. -lldB.IF 

3. -20dB ,IF 


1. -4dB,PLL 

1. -5dB,PLL 

1. -14dB,IF 

560 HZ - 





2. -4dB,PLL 

2. -7dB,IF 

2. -16dB ,IF 

900 HZ 





3. -12dB,IF 

3. -16dB ,IF 

3. -22dB,IF 

- 


Legend : 

1. ms deviation equal to 201 Hz 

2. rms deviation equal to 402 Hz 

3. rms deviation equal to 603 Hz 

IF IF Filter dominant cause of distortion 
PLL PLL stress dominant cause of distortion 
Note: the number indicated is the output SNR in dB relative to the 
SNR for the case in which there is no distortion. 
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Figure 19. PLL Unlock Boundary as a Function of f„ for f 
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frequency (f ) which determines the boundary it is conjectured that 
m 

an approximation to the PLL unlock boundary for the Experimental 
System is given by: 

(Af )(f,) = K(f.) f 2 (4-3) 

rms 2 In 

2 

where Af is the rms deviation, f the PLL natural frequency 
rms n 

(860 Hz) , and K(f^) a function of f^. Equation (4-3) is valid 
only for < .85 f , and is plotted with K(f 1 ) = 1.11 in Figure 19. 

As can be seen there is good agreement between the measured values 
and those given by Equation (4-3). It is beyond the scope of this 
paper to attempt to derive an analytical expression for the unlock 
boundary valid for f 2 > 0.85 f . However, as can be seen in Figure 
19, the unlock boundary seems to approach a constant value of f 2 co « 
In summary then, it is desirable to choose an IF Filter band- 
width and rms deviation so that the output signal undergoes negli- 
gible distortion. The deviation should be small enough so that the 
PLL is not "stressed" excessively. The IF Filter bandwidth should 
be wide enough to attenuate only a negligible percentage of the FM 
spectrum. 

4. 2 Comparison of Theoretical and Experimental Results 

2 

Figure 20 is aplot of 0 as a function of the rms deviation. 

X 

S3 and S4 have essentially the same bandwidth, yet S4 which is 
centered around a higher frequency causes a larger multiplier out- 
put variance. Thus higher frequency components of the signal cause 
2 

larger O and are consequently distorted more than are lower frequency 
components. There is excellent agreement 'between the Experimental 



( volts 
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2 

and Theoretical values (less than 2% error) for G < .25. For 

X 

2 

g > . 25 the agreement is not as good, probably due to the fact 

2 

that the PLL becomes non-linear for larger O . 

X 



Chapter 5 

PLL CHARACTERISTICS IN THE ABSENCE OF MODULATION 

5.1 Introduction 

The purpose of this section is to show that the PLL character- 
istics predicted by the Quasi-Linear Model agree well (within 10%) 
with actual characteristics measured from the Experimental System. 
Since Viterbi f s model is often used (and misused) in determining 
PLL characteristics a comparison is then made between the Quasi- 
Linear Model and Viterbi's Model. It is shown that Viterbi's Model 
is a special case of the Quasi-Linear Model and the regions in the 
PLL parameter space in which each is valid is discussed. 

5.2 Comparison of Theoretical and Experimental Results 

Since there was no junction in the Experimental System from 

2 

which to measure <j>(t), the multiplier output variance (6 ) was mea- 
sured and used to compare to the results predicted by the Quasi- 
Linear Model. Both the theoretical and experimental results are 
plotted in Figures 21 through 23. In each Figure the noise-spectral- 
density (NSD) is held constant while the IF Filter Bandwidth is 
varied. The corresponding PLL signal-to-noise ratio (SNR-^p) f° r 
a wide IF Bandwidth is also indicated as determined from the defi- 
nition [3,4] 

oo 

SNR loop = 10 Io8 10 [A 2 /NSD j |H(s)| 2 fS]. (5-1) 

— OO 

This quantity is commonly used as a measure of the noise in the 
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Figure 22. Multiplier Output and Phase Variance as a Function of IF Filter Bandwidth for No Modulation 
and SNR n = 4.68 dB. 
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VARIANCE/NSD (10 ) 


Theoretical (Quasi-Linear Model) 

— —Quasi-Linear Model with G = A 
X Experimental for 


PLL 3dB bandwidth 




IF FILTER BANDWIDTH = 2g (HZ) 

Figure 23. Multiplier Output and Phase Variance as a Function of IF Filter Bandwidth for No Modulation 
and SNR. =0.2 dB. 
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loop for the case in which the IF Filter Bandwidth is large compared 
to the PLL Bandwidth. It is emphasized that this quantity is included 
just as a matter of convenience and has no meaning when the IF Band- 
width is of the same order of magnitude as the PLL Bandwidth. The 
important parameter is the NSD which is related to the SNR^ oQ p by 
Equation (5-1) . 

In Figure 21 SNI^ = l4.2dB and corresponds to small values 
2 

for the NSD and O' Thus the PLL is operating in the linear region 

<$> 

and replacing G with A in the Quasi-Linear Model gives almost identi- 
cal results. Notice that the error between the experimental and 

theoretical data is less than about 5% even for small IF Bandwidth. 

2 

For small IF Filter Bandwidth the phase variance (o^) is larger than 

2 

multiplier output variance (o ) . As the IF Filter Bandwidth is 

DC 

2 2 

increased, o, approaches a constant value while o increases without 
d> x 

bound. The latter is true since the noise power (variance) increases 
without bound as the IF Filter Bandwidth is increased for a constant 
NSD. Since the PLL does not "track" noise outside of its bandwidth, 
an increase in IF Bandwidth beyond the PLL Bandwidth only adds 
frequency components outside this Bandwidth. Thus the former is 
true since phase error is caused only by frequencies that the PLL 
"tracks." 

The result of increasing the NSD is shown in Figures 22 and 23. 

As the NSD is increased the values predicted by the Quasi-Linear 

Model and the Quasi-Linear Model with G = A differ considerably. 

2 

This is due to the increase in cf^ causing the PLL to operate in a 
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non-linear manner. The Experimental results for O agree well with 

that predicted by the Quasi-Linear Model while differing considerably 

from that predicted by the linear model. This result supports the 

2 

validity of replacing A by Ae 0 as developed in Chapter 1. Figures 

2 

21 through 23 illustrate the complex relationship between o and 

2 2 2 

C^. In general is quite different from.cf^ and in order to deter- 

2 2 

mine <J^ knowing o^, these plots (or the equations used for these 
plots) must be used. 

2 

Figure 24 is a plot of a as a function of CNR for three of 

the IF Filters used in the Experimental System. As can be seen, the 

error between Experimental data and data predicted by the Quasi- 

Linear Model is less than about 20% for all cases. In general, error 

2 

is larger for larger cr and/or smaller IF Filter Bandwidths. This 
is consistent with the assumptions made in developing the Quasi- 
Linear Model (Section 1). 


5.3 Comparison of the Quasi-Linear Model and Viterbi t s Model 

By replacing the equivalent non-linear gain G in Figure 5 by 
A, and allowing the IF bandwidth to become large compared to the 
loop bandwidth, the Quasi-Linear PLL Phase Model becomes similar to 
the model developed by Viterbi [3]. The main difference is that the 
Quasi-Linear Model gives a finite variance for the multiplier out- 
put while Viterbi 1 s model yields an infinite variance (which is of 
course a physical impossibility). The reason for this discrepancy 
is that Viterbi* s model does not account for the fact that in an 
actual system the IF Filter causes the input noise spectrum to have 
a band-pass shape rather than a flat ("white") spectrum. However, 
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the two models do agree on the value for the phase error variance, 

subject to the above mentioned conditions being satisfied. 

Figure 25 illustrates the phase error variance predicted by 

the two models as a function of IF Filter bandwidth for G = A in 

the Quasi-Linear Model. As can be seen, the two models agree to 

within about 2% for BW IF > 10 B pLL « However, the difference becomes 

progressively worse as BW__ approaches zero. For B = 0 the phase 

XF lr 

error variance should be zero since all the noise is filter out by 


the IF Filter. The Quasi-Linear Model prediction coincides with this 

point while Viterbi’s Model does not. As was shown in Section 2.2, 

the Quasi-Linear Model predictions agree well with actual experi- 

2 

mental data. Thus for small (so that G - A in the Quasi-Linear 

Model) and BW TT _ > 10 B Viterbi's model is essentially equivalent 

to the Quasi-Linear Model. The region in the PLL parameter space in 

which each model is valid is illustrated in Figure 26. The boundaries 

2 

illustrated are approximately where the error between predicted 
by the models and actual experimental data is 10%. As can be seen, 
Viterbi’s Model is valid in a subspace of the space in which the 


Quasi-Linear Model is valid. Although Figure 26 was determined for 


no modulation, subsequent data indicates that it is valid whether a 


2 

4 > 


is due to modulation or noise or a combination of the two. 



/NSD (10 ) 



Figure 25* Phase Error Variance Due to Additive Noise of Uniform Spectral Density (NSD). 




Chapter 6 


PLL CHARACTERISTICS IN THE PRESENCE OF 
MODULATION AND NOISE 

6,1 Comparison, of Theoretical and Experimental Results 

Both Experimental and Theoretical Data was obtained for com- 
binations of three different IF Filter Bandwidths, four different 
signals, and three different deviations. The data obtained is plot- 
ted in Figures 27 through 38. In each Figure the IF Filter Bandwidth 
and Signal spectrum is held constant while the deviation is changed. 
The Figure numbers are listed in Table 2 which may be used as an index 
A summary of the maximum errors is also presented in Table 2. The 
first number is the maximum error in dB between the Theoretical 
calculation of output SNR (as determined from the Quasi-Linear Model) 
and the output SNR measured from the Experimental System. The second 

number is the maximum % error between Theoretical and Experimental 
2 

values for 0 . 

x 

As indicated in Section 4.1 the finite IF Filter Bandwidth causes 
output signal distortion for certain combinations of system parameters 
When such distortion exists it is indicated by an asterisk (*) in 
Table 2. Since the effect of IF distortion on the signal was not 
incorporated in the Quasi-Linear Model it is reasonable to compare 
the Model results to Experimental results when IF distortion does 
not exist. As can be seen from Table 2, the maximum error for all 
cases in which IF distortion does not exis.t is 4dB for the output 
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Figure 30. SNR and a x as a Function of CNR for 0 = 5600 Hz, f x * 56 Hz, f 2 » 900 Hz 
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Figure 33. SNR and a as a Function of 
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Figure 36. SNR and as a Function of CNR for 3 = 5600 Hz, f^ = 560 Hz, ^2 ** Hz 
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Table 2 

Comparison of Theoretical (Quasi-Linear Model) 
and Experimental Results 


IF FILTER 




(fl - £2)^21 

5600 HZ 

3300 HZ 

900 HZ 


Figure 27 

Figure 28 

Figure 29 

56 HZ - 

1. ldB,10% 

1. 0. 5dB, 5% 

l.*2dB,15% 

440 HZ 

2, ldB,10% 

2. 0.5dB,10% 

2,*5dB,30% 


3. ldB,10% 

3. 2dB,10% 

3 . *10dB,30% 


Figure 30 

Figure 31 

Figure 32 

56 HZ - 

1. 2dB,10% 

1. 0. 5dB, 5% 

1 4 *2dB, 10% 

900 HZ 

2. 3dB, 10% 

2. 2dB, 15% 

2.*5dB,15% 


3.*12dB,40% 

3**10dB, 30% 

3 . *12dB, 10% 


Figure 33 

Figure 34 

Figure 35 

56 HZ - 

1. 2dB , 5% 

1. 0.5dB, 5% 

l.*5dB,20% 

1250 HZ 

2. 4dB,15% 

2. 3dB,15% 

2 ,*10dB, 20% 


3.*13dB,30% 

3 . *12dB,40% 



Figure 36 

Figure 37 

Figure 38 

560 HZ - 

1. 2dB,10% 

1. 0*5dB,15% 

l.*3dB,20% 

900 HZ 

2. 4dB, 15% 

2**5dB ,40% 

2 ,*5dB, 20% 


3.*10dB,40% 

3.*15dB,40% 



Legend : 

1. rms deviation equal to 201 Hz 

2. rms deviation equal to 402 Hz 

3. rms deviation equal to 603 Hz 

* indicates IF Filter distortion (see Table 2) 

Note: the first number is the maximum difference between theoretical 
and experimental SNR in dB, while' the second number is the 
maximum percent difference between theoretical and 
experimental values of cr^. 
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2 

SNR and 15% for o\ It is emphasized that this is the maximum error 

and that the error is usually much smaller than the above values as 

can be seen from Figures 27 through 38. Table 2 also indicates that 

2 

the error increases as G increases and/or IF Filter Bandwidth 
increases. This is consistent with the assumptions used in develop- 
ing the Quasi-Line ar Model (Chapter 1) . 

6.2 PLL Characteristics as a Function of RMS Deviation 

One of the advantages of using EM is that there is a so called 
m FM improvement" for large modulation indices. Essentially what this 
means is that the output SNR can be improved (increased) by using a 
larger modulation index. This of course assumes that the larger mod- 
ulation index does not cause IF Filter distortion or PLL distortion 
due to loop "stress." Thus there is a "tradeoff" between "EM improve- 
ment" and distortion when the modulation index is increased. 

Figure 27 illustrates the "EM improvement" effect. In this case 
there is no IF or PLL distortion, and an increase in the rms deviation 
(modulation index) causes an increase in output SNR. An increase in 
the deviation by a factor of 2 causes an increase in SNR by a factor 
of 4 (6dB) . Figure 37 illustrates the effect of IF Filter Distortion. 
In this case an increase in rms deviation causes a decrease in the 
output SNR. This is true since IF Filter distortion increases with 
an increase in the rms deviation and for this case the IF distortion 
dominates over any "FM improvement" effect. It should be noted that 
the theoretical values do not conform to the IF Filter distortion 
since this effect was not incorporated into the Quasi-Linear Model. 



From the above discussion it is clear that there is an optimum de- 
viation for a given IF Filter and PLL in the sense that the output 
SNR is maximized. 

6. 3 PLL Characteristics as a Function of IF Filter Bandwidth 
When the IF Filter Bandwidth is of the same order of magnitude 

as the PLL Bandwidth, the noise spectrum is within a frequency range 

that can be tracked by the PLL. Thus the PLL tracks the noise and 

results in a reduced output SNR. As the noise amplitude is increased 

the probability that the carrier will undergo a 2 tt step in phase 

increases [13]. Since the PLL tracks the input noise when such an 

event occurs, there is an impulse noise event at the PLL output. 

This impulse noise event is referred to as a "Type I click" [12]. 

This phenomenon was observed during operation of the Experimental 

System when the* IF Filter Bandwidth was narrow and the CNR was low. 

For wider IF bandwidths, the frequency of occurance of these "Type 

I clicks" decreases which agrees with previous work [12,13]. 

From Figures 27, 28, and 29 it can be seen that a decrease in 

2 

IF Filter Bandwidth results in a larger value for a and a lower 

x 

value for the SNR if the CNR is kept constant. This is consistent 
with the above mentioned fact that the PLL tracks the noise frequency 
components .within the PLL Bandwidth. 

6.4 PLL Characteristics as a Function of the Signal Spectrum 

By comparing the results in Figure 27 with the results in Figure 
36, it can be seen that higher frequency components of the signal 
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undergo more IF and PLL distortion than do lower frequency components. 
This can also be seen by comparing Figures 27, 31, and 33 in which 
the upper cutoff frequency (^2^ increased. Equation (1-19) indi- 
cates that the output NSD is parabolic (i.e. increases in proportion 
to oj ). Therefore, the higher frequency signal components are "imbedded 
in more noise than are lower frequency components. Both the parabolic 
NSD and the IF Filter and PLL distortion effects mentioned above con- 
tribute to the characteristic that in general low frequency signal 
components are demodulated ^ith greater fidelity than are high- 
frequency signal components. 


6.5 PLL Threshold Characteristics 

If the PLL were a linear device the plot of SNR as a function of 
CNR shown in Figures 27 through 38 would be linear. However, due to 
the non-linearity of the PLL the plot departs from that given by the 
linear model for sufficiently small CNR. The point at which the 
actual output SNR is ldB less than the SNR predicted by the linear 
model is termed "threshold" [5]. As can be seen from Figures 27 
through 38 the "threshold" effect is more pronounced for certain con- 
ditions. In general the "threshold" effect is more pronounced (has 
a sharper "bend" in the SNR/ CNR plot) for larger deviations, smaller 
IF Filter Bandwidths, and smaller Signal Bandwidths. 


Quite often an attempt to predict threshold is made by setting 


at equal to some value [1,4,7] and determining the corresponding CNR 
<P 


for the linear model. To check on the validity of this approach, the 
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threshold point was determined for Figures 27 through 38 when the 

"bend" in the SNR/ CNR plot was pronounced enough to give an accurate 

2 

reading. The corresponding values for were then determined from 
the print out generated during runs of the computer program which 
calculated theoretical values. The following observations were made. 


For wide signal bandwidths, threshold occurs at higher values for 
2 2 2 

0^ and a^. The phase variance (a^) at threshold depends mainly on 
the signal spectrum and is fairly insensitive to changes in IF Filter 


Bandwidth and deviation. In, particular, for f^ = 56Hz, and f ^ ” 900Hz 
at threshold - 0.12; while for f- = 560Hz and f 0 = 900H at threshold 

(J) * 1 L Z 

O 2 , “ 0.13. 

9 



Chapter 7 


CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 
7. 1 Conclusions 

The Quasi-Linear PLL Phase Model developed in Chapter 1 is 
realistic and useful in the sense that it accurately predicts phase 
error and multiplier output variances as well as output SNR of an 
actual PLL. It is an improvement over previous models in that it 
accounts for (in a statistical sense) the fact that the IF Filter 
Bandwidth may be narrow compared to the PLL Bandwidth and that the 
phase error may be large enough to cause operation in a non-linear 
region. The use of a realistic modulating signal contributes to 
the usefulness of the results of this paper. For the second-order 
PLL and band-pass type modulating signals used in the Experimental 
System the error between actual and predicted data was less than 
15%. However, there is no reason to believe that a similar bound 
on errors would not be found for different order PLL’s and modu- 
lating signals. 

The following observations were made concerning the occurence 
of threshold. The threshold effect is more pronounced (has a sharpe 
"bend” in the SNR/ CNR plot) for narrow IF Filter bandwidths, for 
larger modulation indices, and for narrow signal bandwidths. The 
phase error variance at threshold depends mainly on the modulating 
signal spectrum, and varies considerably as the signal bandwidth or 
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center frequency of the signal spectrum is changed. The above 
threshold information is important since it is desirable to operate 
the PLL above threshold. 

7. 2 Recommendations for Future Work 

One of the difficulties encountered in using the Experimental 
System was that the VCO "phase jitter" made measurements of the out- 
put SNR above approximately 18dB inaccurate. The use of a voltage- 
controlled crystal oscillator might reduce this inaccuracy and allow 
measurements at larger SNR's, and variances. Then threshold character- 
istics for large modulation indices could be investigated. 

An important effect encountered when making measurements on 
the Experimental System was output signal distortion caused by the 
finite IF Filter Bandwidth. It would be very helpful if an analytical 
expression could be developed which gives the relationship between 
IF Filter distortion as a function of the IF Filter Bandwidth, modu- 
lation index, and signal spectrum. 

The theoretical and experimental techniques developed in this 
paper were applied to a second-order PLL. These same techniques could 
be applied to PLL*s that are other than second-order. One particular 
PLL which has received attention [10] to which these techniques may 
be applied is the Multi filter PLL. 
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Appendix I 


SPECTRAL DENSITIES OF THE PLL 


Let ©i(s), 02 (s) » $(s), N(s), Y(s), and O(s) represent 
the Fourier Transforms of Oi(t), © 2 (t), $(t)» n(t), y(t), and 
o(t) respectively. Then, from Figure 5, 

«(s) = ©l(s) - 02 (s) 

= 01 (s) - [N (s) + G *(s)]. 


or $(s) J 
Thus $(s) = 


Msi + q 

S© 1 (s) 
s + GK F(s) 


= 01 (s) - K ^ (s) N(s). 

K F(s) N(s) 
s + GK F(s) * 


Since H(s) 


GK F(s) 
s + GK F(s) 


and 1 - H(s) 


s 

s + GKF(s) 


(x-i) 


(X-2) 


Equation 1-1 becomes 

i(s) = [1 - H(s) ] 01 (s) - — H(s) N(s). (1-3) 

Since i(t) and n(t) are independent 


$(io) = |l - H(s) | 2 0i(u) + | |H(s)| 2 N(u). 

G ' 


Again from Figure 5 

Y(s) = N(s) + G $(s). 

substituting for 0(s) as given by Equation (1-3) 
Y(s) = N (s) + G[{l - H(s) }©i (s) - H(s) N(s>] 

- G[1 - H(s) ] ©i (s) + [1 - H(s) ] N(s) . 


(1-4) 


(1-5) 


( 1 - 6 ) 



Thus the spectral density is 
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Appendix II 


BAND-PASS FILTER TRANSFER FUNCTION 


Consider first the low-pass section of the Filter. Since 
Equation (2-1) is the transfer function squared it is necessary 
to find the polynominal V^(s) such that 


64e‘ 


V s > V 4 ( - s) 


2 2 

1 + e T^(s) 


(II-l) 


64e‘ 


1 + e 2 * [64s 8 - 128s 6 + 80s 4 - 16s 2 + 1] 


1 + e 2 T^(s) = 0. 


Thus one must find the roots of 

(II-2) 

The procedure to do so is given by [6] and yields the following 
4 roots in the left-hand complex plane 


r i = - a i + jb i 
r 2 = " a l " - 1b l 
r 3 = -a 2 + j fc 2 
r 4 = -a 2 ” Jb 2 


where a^ = 0.203721927 
b ± = 1.046636712 
a 2 = 0.491828394 
b 2 « 0.433531121 


(II-3) 


Therefore 


V 4 (s) = (s +a 1 -jb 1 )(s +a x +jb 1 >(s +a 2 -jb 2 ) (s +a 2 +jb 2 ) 


= s 4 + 2(a^ + a 2^ s ^ "** "** ^2 "** 4a^a 2 )s 2 


+ 2(a,A 0 + a 0 A n )s + A^A 2 , 


12 


2 1 


(II-4) 


where 


2 2 2 2 

A^ = a^ + b^ , and A 2 = a 2 + ^2 
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Substitution of Equation (II-3) into Equation (II-4) yields 
V 4 (s) = s 4 + f 3 s 3 + f 2 s 2 + f 1 s + f 

where 


f Q = 0.4887118438 


f 1 = 1.29350660 


f 2 => 1.96758027 


f 3 = 1.39100322 


(II-5) 


Thus the transfer function for the low-pass section is 

1 


B(s) 


LP 


(k) +f 4J +f 4 2 ) +f 4 2 ) 


+ f 


0 


(II-6) 


where tt 2 is the upper break frequency of B(s). Likewise the transfer 
function of the high-pass section is 

1 


B(s) 


HP 


(?) * ■£) * ',(?) * 


+ f 


0 


where is the lower break frequency of B(s). 


(II-7) 


The overall transfer function (B(s)) is 
B(s) * B (s) Lp Bts)^ 


8 . 7 , 

g g s + g ? s + 


(II-8) 


.+ g 2 s +g x s + g 0 


where 


>0 ~ f o w i 


, * 

- _L f f 3 

8 1 " + f 0 f 3 W l 



+ 
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